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The non-modal kinetic theory of the kinetic drift instability of plasma shear flows 
[Phys. Plasmas, 18, 062103 (2011)] is extended to the investigation of the long-time 
evolution of the hydrodynamic ion temperature gradient and resistive drift instabil- 
ities in plasma shear flow. We find, that these hydrodynamic instabilities passed in 
their temporal evolution in shear flow through the kinetic stage of the evolution. In 
linear theory, this evolution involves the time dependent, due to flow shear, effects of 
the finite Larmor radius, which resulted in the non-modal effect of the decrease with 
time the frequencies and the growth rates of the instabilities. 
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I. INTRODUCTION 



The investigation of the long-time mode of behaviour of the regimes of tokamak plasma 
improved confinement, where E x B flow shear stabilization mechanisms^ plays a key role, 
is the theme of major importance in the physics of continuous tokamaks operation. Gy- 
rokinetic theory- gives a foundation for the investigating microinstabilities, which cause the 
anomalous transport in fusion plasmas. While there has been significant progress in gyroki- 
netic treatment of tokamak turbulence (see, for example, Refs.-^-^), there key issues still 
remain in application of the gyrokinetic theory to analytical investigations of the long-time 
evolution of plasma shear flows turbulence. 

In Ref.— we have developed new kinetic description of the plasma shear flows, using the 
Kelvin's method of shearing modes or so-called non-modal approach as its foundation. The 
governing equation in that theory is the integral equation for the perturbed electrostatic 
potential, which is the formal solution of the initial value problem for Vlasov-Poisson (V-P) 
system. The important procedure in its derivation is the transformation of V-P system to 
sheared (in space and velocity) coordinates convected with flow and accounting for by this 
mean the effect of spatial time-dependent distortion of plasma disturbances by shear flows. 
In that equation, velocity shear reveals as the non-modal time-dependent effect of the finite 
Larmor radius. Because of the shearing of perturbations in shear flow, the component of 
the wave number along the direction of the velocity shear experiences the growth with time, 
and therefore results obtained for the hydrodynamic drift instabilities on the base of the 
fluid description are valid only for finite times till kj_pi <C 1 is valid ( pi is thermal ion 
Larmor radius). Therefore, not only kinetic, but also typically hydrodynamic instabilities 
of shearless plasma require the kinetic description for plasma shear flows for proper treating 
their evolution on times, at which k±pi becomes not small and approaches unity and above. 
It reveals^, that conventional kinetic theory, as well as the gyrokinetic theory, in which 
the time dependence of the perturbed distribution functions and fields are considered in a 
canonical modal form, ip ^ exp (—iut), are valid only at the initial stage of the shear flow 
turbulence evolution, at times less then the inverse velocity shearing rate (Vq) -1 . The theory 
developed displays that the experimentally observed suppression of drift-type turbulence and 
the improved energy confinement are appreciable at t > (V^)" 1 > (7)" , i.e. just after that 
time. In our analysis of the kinetic drift instability we have obtained^, that in linear theory 
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of the kinetic (universal) drift instability, shear flow leads at time t > (Vq) _1 to the non- 
modal decrease with time the frequency and growth rate and to ultimate suppression of that 
instability 

The basic transformation of the V-P system to sheared coordinates, which resulted in 
the removing from the Vlasov equation the term, which contains the spatial inhomogeneity 
introduced by the velocity shear, is presented and discussed in Sec. II. In Sec. Ill, we 
apply the linear non-modal kinetic theory of Ref.—, extended onto the accounting for the 
inhomogeneity of the ion temperature, to the investigation of the temporal evolution of the 
hydrodynamic ion temperature gradient drift instability of plasma shear flow. We find that 
the ordinary modal theory of this instability is valid only for times less than the inverse 
velocity shearing rate. At times t s 3> t > where t s = (V^kyPi)' 1 , the non-modal 

effects of the decrease of the frequency and growth rate develop. At the final stage, at times 
t t s the perturbations of the cell type with zero frequency occurs. 

In Sec. IV, using the integral equation for electrostatic potential obtained in Ref.—, ex- 
tended onto the accounting for the collision of electrons with neutrals, we consider temporal 
evolution of the resistive drift instability of plasma with comparable ion and electron tem- 
peratures, as it is at the edge layer of tokamaks. We came to the same conclusion about 
the importance of the kinetic theory for the proper treating of the resistive drift instability 
in shear flow. It is shown, that the reducing with time the frequency and growth rate is a 
common future of the evolution of the electrostatic drift instabilities in plasma shear flows. 
A summary of the work is given in Conclusions, Section V. 

II. VLASOV-POISSON SYSTEM OF EQUATIONS 
IN SHEARED COORDINATES 

We start with the Vlasov equation for species a (a = i for ions and a = e for electrons), 
immersed in crossed spatially inhomogeneous electric field, Eo (r) and homogeneous magnetic 
field B || e z , 




We use a slab geometry with the mapping (r,6,ip) — > (x,y,z) where r,9,<p are the radial, 
poloidal and toroidal directions, respectively, of the toroidal coordinate system. In this 
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paper, we consider the case of plasma shear flow in linearly changing electric field, Eo (r) = 
(dEo/dx) xe x with dE^/dx = const . In that case 

c OE 

V (r) = V (x) e y = -^^ xe y = V o% e y ( 2 ) 

with spatially homogeneous, Vq = const, velocity shear. The possible spatially homogeneous 
part of shear flow velocity is eliminated from the problem by a simple Galilean transforma- 
tion. It was obtained in Ref.— , that transition in the Vlasov equation from velocity v and 
coordinates x, y, z to convected coordinates v in velocity space, determined by 

v x = v x , v y = v y + VqX, v z = v z (3) 

and to sheared with flow coordinates x, y, z in the configurational space, determined by 

x = x, y = y + V^tx, z = z (4) 

(it is assumed that inhomogeneous electric field, and respectively shear flow originate at 
time t = i(o) = 0) transforms the linearized Vlasov equation for the perturbation of the 
distribution function f a = F a — F 0a , with known equilibrium distribution F 0a , to the form, 

9fa , df a , , T// .. df a , df a , 8f a , , df a 



+ V ax-^ + {Vay ~ V ax V' a t) — - + V az —^ + U ca V ay —^- - (u ca + V£) V ax - 



dt "* dx v y a ' dy w dz a y dv ax v a ' "-dv. 



ay 



e a f dip y/ t di f \ dF oa ! e a dp dF 0a | e a dp dF 0a 
m a \dx a dy ) dv ax m a dy dv ay rn a dz a dv az 

(u c is the cyclorton frequency of ion (electron)) which is free from the spatial inhomo- 

geneities originated from shear flow (see also Eq.(8) in Ref.— ). The Fourier transformation 

of Eq.(j2]) over spatial coordinates x,y,z with the electrostatic potential <p(r,t) determined 

function of coordinates 

P (x, y,z,t) = J p (k x , k y , k z , t) e^ x+lk yy + ^ z dk x dk y dk z , (6) 

gives the equation for the separate spatial Fourier harmonic of f a with wave numbers 
k x , ky, k z 

() ' ,a + (i (k x - V^tky) v ax + ikyV ay + ik z v az ) f a (v Q , k, t) + u ca v ay ^- - (u ca + V£) v ' ln 



i \" \'^' ' u y / ax i v,u y ^ ay i " ,tj z ^ olz / j oc \ * en "? " j 1 w cot ^ ay \ w cot 1 'a/ ua. 



i—tp(k,t) 



\K X VnlKy) — h %Ky— r 1rC z ~pz 

dv ax dv ay dv az 



(7) 



in which wave numbers k x ,k y ,k z and time change independently. It follows from Eqs.flS]) 
and (J7|), that the transformation of the Vlasov equation to convected-sheared coordinates 



([3]), OH) converts the spatial inhomogeneity into the time inhomogeneity. That prevents the 
application of the spectral transforms in time to Eq. (J7J) and to obtain the ordinary dispersion 
equations valid for any times. It reveals^ that the ordinary modal solution to Eq.(j7]), for 
which ordinary dispersion equation may be obtained, exists only at time t <C (^o') _1 ; solution 
becomes non-modal^ with time dependent frequency and growth rate for larger time. With 
coordinates x, y, z of the laboratory set of references, transformation has a form 

p(x,y,z,t)= f p(i,,^t,t|t iMt, ''(w' s )* i JWl,4, 



<p (k x , k y , k z , t) ei k *- v ° tk «) &+ik ^ +ik **dk x dk y dk z . (8) 

It follows from Eq.flS]), that separate spatial Fourier mode in convected-sheared coordinates 
becomes a sheared mode with time dependent wave number k x — V^tky in the laboratory 
frame. 

Usually, however, only the transformation to the convected coordinates ([3]) in velocity 
space, without the transformation to the sheared coordinates (j4]) in the configuration space, 
is used in the kinetic theory of plasma shear flows^ 1 ^. After such transformation the 
linearized Vlasov equation becomes 

df a ,,df a ^df a df a . „ 8f a 

Get /a ,\ dF() a , , 

= 5£ v "<'■*> "e»- (9) 

Traditionally^ 1 ^, in proceeding with derivation of the governing equation for f a , the spatial 
Fourier transform in the laboratory configuration space, 



for the electrostatic potential (p and for f a is adapted with assumption of the "slow variation" 
of Vo (x) with spatial coordinates. That gives the following equation for f a \ t, k, vj : 

df a ( * \ df a df a 

ik y V {x) fa-i (kir) fa+UJ ca VayT^ ~ {u) ca + Va)Va X - 

dF ( 



fu "Tur/Jii l J J ix • ~ca~ay p. \~ cu . ■ a/ a 

UL \ / UUax '-"Jay 



= -^(k,t)^. (11) 

In Eq. ffTO]) . the flow velocity shear reveals in the formation of elliptical orbits of particles, 
with velocity coordinates v±, <p 15 

v x = v ± cos0, v y = yfrjv± sin0, <p = (pi - y/rju c t, v z = v z , (12) 



where rj = 1 — Vq/u c , and with modified gyro-frequency y/rju c . That effect, however is 
negligible for \Vq\ <C 0J C i- Then, velocity shear is absorbed into the identical for both plasma 
species Doppler shifted frequency^^ Cj = u — k y Vo (£), and, in fact, is excluded from the 
subsequent analysis. 

Let us analyse the results of the application of the Fourier transformation ( flOl) to Eq.flH]) 
without the assumption of the "slow" spatial variation of the flow velocity. That gives the 
following differential equation in wave-number space: 

~dr~ V ° ky K~ l ( > fa UcaVay d^ " {Uca + K) Vax d^ y 

= -^(k,t)^. (13) 
m a V / (TV 

It is interesting to know the relation between solutions of Eqs. (j7]) and fflBl for f a . For 
the receiving from Eq. tfTBl the equation which couples f a and (p of the separate spatial 
Fourier mode, as it is in Eq.(jSJ), we have to exclude from Eq. fflBl the differential operator 
—Vnky^-, due to which the Fourier mode of f a appears to be coupled with all Fourier modes 
of the electrostatic potential and depends on the integral of ip over wave-number space. The 
characteristic equation 

dt = -*h- (14) 

V'k 

gives the solution k x + V^tky = K x , where K x as the integral of Eq. fflB"]) is time independent. 
It reveals that f a = f a \^K x ,k y ,k z ,tj = f a \k x + VQtky,k y ,k z ,tj, i.e. the wave number 
components k x and k y have to be changed in such a way that k x + V^tky leaves unchanged 
with time. The solution to Eq. fl!3j) for f a can't be presented in the laboratory coordinates 
in a form, in which the time and spatial dependences are separable, as it is for the normal 
mode solutions of Eq. fflTT) obtained with assumption of the "slow" spatial variation of 
the flow velocity^. If we use, however, k x = K x — V^tky in Eqs. flTU]) and (113"]) . we obtain 
for the electrostatic potential the presentation (jSJ), and we obtain Eq.([7j) for f a , with time 
independent K x = k x , k y = k y , k z = k z . The obtained results prove, that the solution of the 
Vlasov equation in the form of the separate Fourier harmonic with time independent wave 
numbers may be obtained only in convected-sheared coordinates. That solution reveals in 
the laboratory frame as a shearing mode jSJ) with time dependent x-component of the wave 
number. 
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So, we have two procedures for the proper performing of the spatial Fourier transformation 
of the Vlasov equation for plasma shear flow, which give the same result. The first one is 
to apply at first the transformation §2§ to convected-sheared coordinates, and then, to 
perform the Fourier transform of the Vlasov equation over spatial coordinates with time 
independent wave numbers k x , k y , k z . The second procedure is to come to these time 
independent wave numbers through the solution of the characteristic equation Eq. fH]) . when 
the transformation to sheared coordinates in configuration space in the Vlasov equation does 
not perform. The sheared flow leads to the observed in the laboratory frame continuous 
inclining with time of the plane of constant phase the waves. That reveals in the observed in 
the laboratory frame the time dependence of the wave number k x = K x —V§k y t in the solution 
for the f a . The oversimplification of the problem, which resulted from the application of the 
assumption of "slow spatial variation of Vq (x)" , leads to the overlooking of that principal 
effect of shear flow. It is obvious, that the time dependence in K x may be neglected only in 
the case of negligible velocity shear, or when the very short evolutionary time is considered. 
Really, for k y ~ k x and V ' ~ 7 for time t > 7 _1 the we have kyV^t > k x in the integral 
K x = k x + V^kyt. Therefore the assumption of "slow spatial variation of flow velocity" is not 
valid for the investigations of the effects of shear flow in real experiments, where observed 
velocity shearing rate may be of the order or above of the growth rate of the instability and 
the time of the observations is of the order of the inverse growth rate or longer. 

With leading center coordinates X, Y, determined in convective- shearing coordinates by 
the relations 

x = X — sin 0, y = Y H — cos <p + Vnt (X — x) , 

y/rjujc rjujc 

z\ = z- v z t, (15) 

the Vlasov equation (JTJ, in which species index is suppressed, transforms into the form 13 

dF ^ e ( dcp dF dip dF 



dt m^/fjLO c \dX dY dY dX 
| e ^/rjuJc ( dp) dF dp> dF \ e dip dF _ Q 
m v± \ d(pi dv± dv± d<p\ ) m dzi dv z 

in which any time dependent coefficients are absent. It follows from Eqs. ffTo"]) . that in shear- 
ing coordinates a particle gyro-motion is different from the ones in convective coordinates, 
determined by Eqs. (|T2|) . Now it consists in the rotation with modified gyro- frequency and 
stretching of gyro-orbit along coordinate y with velocity —Vq (x — X), which is negative for 



x > X and is positive for x < X. With leading center coordinates (fl5|) the Fourier transform 
(HO) becomes 



ip (x, y,z,t) = J (f (k x , k y , k z , t) exp 
v± 



ikr I X — 



sm i 



+ik y I Y + 



rjoj c 



cos0 I +tV't^±~ : 



sin d> + zA^z 



dk^dkydk % 



J ip (k x , k y , k z , t) exp [ik x X + ik y Y + ik z 



. k y v± COS . — V^tA: y ) V± COS0 



dk^-dkydk % 



. fei (*) v± . 



sin (0! - y/rju ci t - 9 (t)) 



dkr^dkydkg 



where 



k\ (t) = (k X - V^tkyf + ~k\ 



(17) 



(18) 



and tan# = k y j \frj{k x — V^tky). It follows from Eq. ffTT|) that finite Larmor radius effect of 
the interaction of the perturbation with time independent wave numbers k x , k y , k z with ion, 
Larmor orbit of which is observed in sheared coordinates as a spiral continuously stretched 
with time, appears identical analytically to the interaction of the perturbation with wave 
numbers kx — V^tky, k y j ^/rj, k z with ion, which rotates on the elliptical orbit that is observed 
in the laboratory frame. The time dependence of the finite Larmor radius effect is the basic 
linear mechanism of the action of the velocity shear on waves and instabilities in plasma 
shear flow. 

In what follows, we consider the equilibrium distribution function F i0 as a Maxwellian, 

(19) 



n (X) ( v 2 + v 2 



(2tto 2 ^ 3/2 



exp 



assuming the inhomogeneity of the density and temperature of plasma shear flow species 
on coordinate X. In this paper we assume, that velocity shearing rate Vq is much less than 
the ion cyclotron frequency u ci , and put rj = 1. The solution of the Vlasov equation for 
the perturbation / (£, k x , k y , k z ,v±, (ft, v z ,zi) of the distribution function F, f = F — F with 
known F$ is 



s= e - 

m 



1 dip dF uj c dip dF dip 8F 



uj c dY dX v± d(pi dv± dz\ dv z 



dt'. 



(20) 



to 



Using solution f l20|) for all plasma species in Poisson equation for the potential tp (r,t), 

A <p (r, t) = -4tt e a J fa (v, r, t) dv a , (21) 

a=i,e 

we obtain integral equation^ 1 ^, which governs the temporal evolution of the separate spatial 
Fourier harmonic of the electrostatic potential <p (k, t) in plasma shear flow and is capable 
of handling linear as well as nonlinear— evolution of the electrostatic instabilities of plasma 
shear flows. 



III. HYDRODYNAMIC ION TEMPERATURE GRADIENT INSTABILITY 

In this section, we consider the temporal evolution of the hydrodynamic ion temperature 
gradient instability in plasma shear flow. By using the methodology stated in Sec. II, we 
obtain the integral equation for the electrostatic potential $ (k, t) = ip (k, t)Q(t — to), where 
© (t — to) is the unit-step Heaviside function (it is equal to zero for t < t and equal to unity 
for t > t ) 

t 



J dh-^- {$ (k, tj) [- (1 + T) + A 0l (t, h)}} 

to 

t 

j dt 1 ^-{^(k,t 1 )A 0l (Mi) (l-e-**^*-*) 1 



So 



t 

+ J dt&fctJAK (Mi)e-^« (t - tl)2 



to 



X 



ikyV d i - iu T i - k 2 z VTi [t - tij - -u T ik 2 z VTi (t - h)' 



t 



+ iu Ti J dt^ (k, h) A u (t, h) e'h&hit-tif 



to 

t 



+ T jdt, + ik y v^ (k, to) e-h^t-t,? = 0> (22) 

to 

where T = T»/T e , UTi = kyVdid In T^/d In n i? t> di = (cTi/eB) dlnriio/dx is the ion diamagnetic 
velocity, 

A 0l (t, h)=Io (k ± (t) k ± (h) pi) e - W(^W4*L(*i)) ; (23 ) 



t 








T 



S 



FIG. 1. The sequence of the characteristic times for Eq. (|22p for the long wavelength perturbations 
with k± (to) pi < 1. 



Eq. ( 122ft is the extension of the Eq.(25) in Ref.-^ onto the accounting for the inhomogeneity 
of the ion temperature. Naturally, it is not possible to obtain explicit single analytical 
presentation of the solution to Eq. (!22|) . which is valid for any desired time. The exceptional 
advantage of the application of the non-modal approach, which uses the wavenumber-time 
variables, is a possibility to perform the analysis of the electrostatic potential evolution at 
any finite time domain and with an arbitrary initial time to- 111 our analysis, initiated in^, 
we distinguish different characteristic times of the evolution process of the long wavelength 
perturbations with k± (to) Pi <C 1: to is the time of the perturbation origin; t( ) = is the 
time of the electric field (shear flow) origin; times t(i) = (Vq') -1 , t( 2 ) — t s — (V^kyPi)^ 1 define 
different stages of the non-modal evolution of the potential under the action of the shear 
flow. In these time intervals different tasks in the solving of Eq.( )22|) arise (as it was with 
Eq.(25) in Ref.^). We will derive the linear solution of Eq.fl52J) at time (VJ) -1 > t > t of 
the modal instability development, at which shear flow effects are underdeveloped. Then, 
we will derive the solution at time t ^> (Vq')~ , at which the non-modal effects of the shear 
flow become pronounced. And, at last, we will derive solution at the final stage of the 
perturbation evolution at time t ^> t s — (V^kyPi) 1 ^> (Vq') -1 , at which even initially long 
wavelength perturbations with k±pi < 1 at time f<t S) becomes the short wavelength ones 
with k±pi 3> 1 at t 3> t s . 

The long-wavelength drift perturbations with kj_pi < 1 at time t = 0, will be the long 
wavelength ones at time t < t s = (V^kyPi)' 1 . For this limit, we use the approximations, 



x 




) 



(24) 
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k + P \ ( hkX (t + h) - h 2 y (v f (t 2 + t 2 ) ) o < / , } 



(25) 



A u (t,ti) wl-p? (*1 (t) + kl (h)^ 
b a + p 2 (2k x k y V£ (t + *0 - k 2 y (V^f (t 2 + 1 2 ) ) (h) 



(26) 



where b{ = 1 — A: 2 p 2 , &a = 1 — 2k 2 L p 2 , k]_ — k 2 + k 2 , and 9 (t) indicates that the shear flow 
emerges at t — 0. In Eq. (122 jl . we use the approximations 1 — exp [—{l/2)k 2 v 2 ^ i (t — ti) 2 ) — 
—(l/2)k 2 Vj< i (t — ti) 2 , that corresponds to the weak ion Landau damping. Assuming that 
the magnitude of the initial potential $ (k, t — to) is negligibly small, the integral 



to 



t 



dt x ^U (MO (l-e-^>Ti^"' 2 



-fc 2 4, (t-to) 2 $ (k,t 



(27) 



may be ignored in Eq. (l22|) . This approximation is reasonable only for the exponen- 
tially growing, faster than (t — to) or (t — to) 2 , potential $ (k, t); then $ (k, t ) (t — t ) 
and (<i$ (k, t) /dt) | t=to • (t — t ) 2 are exponentially small with respect to $ (k, t). We assume 
that k 2 z v\ e (t — ti) 2 ^> 1, that corresponds to adiabatic electrons, and neglect exponentially 
small electron terms in Eq. (I22[) . In result, we obtain the simplified integral equation, which 
describes the linear evolution of the electrostatic potential at time t <^t s , 



(T + k\p 2 ) J dt 



to 



d 3 ^(M0 

dtf 



to 



, rf 2 ^(k,tl) r, , , ,9 9/ 

~d~t\ L * ' yVdi ~ zVTi ^ x ' 



iuj T ik\p 2 



r rf 2 *(M0 
dt 2 



p 2 / d*i 



dti 



Kk y v ' (t + to - ife 2 (v ') 2 (t 2 + 1\) ) e (to 
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+ Pt J dti 

to 

x 



(28) 



where new variable determined as d 2 ^/dt 2 = $, is introduced. Now we obtain the 
solution to Eg. (128]) in the limit of large value of the parameter rji ^> 1, at which long 
wavelength with k±pi < 1 hydrodynamic ion temperature gradient instability is developed 
in a shearless plasma. By integration in parts of Eq.(|28l) and neglecting the initial values of 
dty (k, t) I dt and \& (k, t) at t — to with assumption of the exponential growth for \l/ (k, t) at 
time t 3> to, w e obtain in the zero approximation the equation 
t 



dti 



, , d 3 ^ (k, ti) , 
(T + k 2 ± p 2 ) + iu Ti k 2 z v Ti b^ (k, ti) 



0. 



(29) 



to 



in which the dominant terms in the left-hand side of Eq.( l28|) are retained, and the right- 
hand side, proportional to k 2 p 2 <C 1, is neglected. The solution to Eq.( l29~]) is ^/ (k, t) = 
C exp (— iu (k) t), where the frequency uj (k) is determined by the known equation 



u 3 (k) = -u Ti kyJ l klpf) 



(30) 



(t + m) ■ 

The accounting for the small right-hand side of Eq. (128]) will modify that solution. We seek 
in the next approximation the solution for \l/ (k, t) in the form 



#(k,t) = Cexp(-zw (k)t + er(k,t)) 



(31) 



where the higher order correction, <r(k, t), we find by use the procedure of the successive 
approximations 13 , which gives the following equation for a (k, t): 



t 

3 (r + k 2 ±P 2 ) u 2 (k) J dM?0Mi)^ 



to 



-p! J dt x ± [ [k x k y v> (t + to - \k 2 y (v*f (t 2 + e (to 



to 



i -^Ti 2 / i+ 

% ~ Pi J 1 

to 



d ^^ tl] (Wo' (t + to - \k 2 u (v>f (t 2 + tl) ) e (to 

xk 2 z v 2 Tl (t-tO 2 , 



(32) 
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t t t t 



FIG. 2. The domain of the integration [to,t] over time in Eq. (|36p for the perturbations with 
k± (t ) Pi » 1. 

where ex (k, t) is neglected in ^ (k, ti) in the right-hand side. The solution to Eq.( l32l is 
obtained straightforwardly and is equal to 

a (k, t) = -h^X [V^tf + -klp\u (k) tX (V^t) 2 , (33) 

where 

1 - T - 2k\p\ 



X 



(i - kip?) (t + kip*y 

This result displays, that linear effect of the shear flow at time t s ^> t ^> (Vq)~ consists in 
the non-modal decrease of the frequency and the growth rate with time. Qualitatively the 
same result was obtained for the drift kinetic instability in Ref.— . 

It may be anticipated that the most substantial effect of the instability suppression will 
be attained at time t approaching t s and that suppression will continue at t > t s . For that 
time, however, solution to Eq.( j22l) may be obtained only numerically. For time t 3> t s , we 
have k± (t) pi ~ kyV^tpi — t/t s 3> 1 and small parameter t s /t ~ t s /to 1 will appear in 
Eq. (l22l) . that admits the receiving of the analytical solution^ to Eq. (l22l) . For that time 



A ^~7m e * v -*•<■*-**)■ (35) 

For the application the asymptotics fl34|) . f )35|) for time t\ in the whole interval t > t\> to, 
we have to consider in integral equation (13"!?|) the initial time to > t s . At that case the 
condition (p (k, t) = 0, for (V'q) -1 < t which was used for the function $ (k, t) in Eq.f l22p . in 
which time to was considered as the time preceding the development of the modal instability 
and appearance of the shear flow, is not applicable now, and we have to restore the equation 
for (p (k, t), obtained in Ref.=. It is important to note also, that for the decaying with time 
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solution for (p(k,t), the approximation of the neglecting by initial value (k, t ) used in 
Eq. fl28l) is not justified and we have to retain these terms in equation considered. We obtain 
finally the following integral equation: 

t t 



;i + t) / dtt 



dip (k, t) 



dt 



t 



dti 



dt t 



<p(k,t) 



to 



-.e 2 



^(^-^l) 2 



+ / dtnp(k,t) 



to 



-.e 2 



(*-*i) 2 



x ( ikyV d i - -uJTi^i (t - h) - k z v Ti (t - ti) 



+ V 9 (Mo) 



t V27T 



Mt-to? 



- 1 



(36) 



where = t s 2 + k 2 z v\ { . At t S> t s , the right-hand side of Eq.( 13"6"j) is proportional to 

small parameter t s /t ~ t s /to. In zero order in that parameter we have the equation 

t 

f dti(d<f(k,ti)/dti) = with solution <p(k,t) = const. The solution, which accounted 

to 

for the terms of the order of t s /t in Eq. (l36|) . we find in the form (f (k, t) = C exp (a (k, t)), 
where a (k, t) = O (t s /t). Omitting a (k, t) in the right-hand part of Eq. (l3"61) . we obtain the 
approximate solution for tp (k, t), 



Lf (k, t) = (fo exp 



kyVdl 

t(i+T) y 2k { 



2 klvl; 



2 
7T 



(37) 



which is qualitatively the same as it was obtained in Ref.— for the kinetic drift instability 
for large time, t 3> t s and is far from ordinary modal solution. 



IV.NON-MODAL KINETIC EVOLUTION 

OF THE HYDRODYNAMIC RESISTIVE DRIFT INSTABILITY 

In this section, we consider the temporal evolution of the resistive drift instability in 
plasma with comparable temperatures of ions and electrons. In convective electron guiding 
center coordinates, linear perturbation f e of the steady state electron distribution function 
F 0e is governed by linearized Vlasov equation with number-conserving Bhatnagar-Gross- 
Krook (BGK) collision terrn^ 

dfe _ e dip dF 0e e u ce dip dF 0e ^ e dip dF 0e _ y _ n ei (r,t) 

dt m e u ce dY e dX e m e v± d<fii dv± m e dz dv z e \ e n e0 
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In Eq.( l38|) n e i is the perturbed electron density, and X e is electron leading center coordinate. 
The solution to Eq. fl38|) is 

t 



fe 



1 dip dF e0 u ce d(p dF e0 dip 8F e0 



to 



u ce dY e dX e v± d<fii dv± dz\ dv 
t 

+ u/[e-»™*±^F 0e dt 1 . 
J n e0 



(39) 



Using solution ( 13"9~|) with Maxwellian distribution F 0e , f[T9~j) . and solution (}2"0"]) for the pertur- 
bation of the ion distribution function in Poisson equation (|2T]) . we obtain in the quasineu- 
trality approximation the following equation for the potential $ (k, t) = <p (k, t) 6 (t — t ) 
for low frequency, dip/dt <C u C iip, perturbations of drift type: 

t 



dtx <i (1 + T) (k, to - J- ($ (k, to ^ (f , *0) e 



-ifcX.(t-ti) 2 



to 



z / dti$ (k, tO /c^A 0i (t, tO + TG e (k, t) 



(40) 



where T = Tj/T e and O (t — to) is the unit-step Heaviside function. G e (k, t) determines the 
nonadiabatic part of the electron density perturbation, 



™ e i(k,t) = -^$(k,t) + ^G e (k,t). 



(41) 



and is equal to 



t 

G e (k, t) = J dt x (Jj- ($ (k, tO e-"-**-*)) + ikyVte* (k, tO e 



to 



l2 / ,7/ ., ,1, -/- > „-f e (t-ti)-§fef4 e (t-ti) 2 



(42) 



+ 47rz/ e A De y dtiriei (k, tO e" 

to 

Now we apply the methodology developed in Ref.— to calculate the approximate solution 
to system (HO"]) . fj4*2|) for long wavelength perturbations with k± (t) < 1 for times limited 
by the condition (V^') -1 < t < t s . Making use the approximation (123]) we present Eq. (j4T)]) in 
the form 

t 



to 
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f , /d$(k,ti) ., ^ „ A 

r + ^/ <iti (^r^ +! ^'* (k -' i) Jv "V 2 

to 



/ i ./. .. /- i '1-exp k 2 z v 2 Ti {t-t x f 



d$(k,h) , .. _ ,_ A fk x (t + h) (t 2 + tj 



i y , /d$(k,ti) ., 







+ rT*M 7 ** (k - h) wM~ Vit > + tTW; G < (k ' () ' (43) 



where u (k) is 

^ ( ^ = l + fc^ Vde. ( 44 ) 

p s is the ion thermal Larmor radius with electron temperature, and Vde is the electron 
diamagnetic velocity. The first term in the right-hand side of Eq. (l43j) determines the ion 
Landau damping; this term is the same as in plasma without shear flow. The next three 
terms originate from shear flow and determine the corrections to the frequency and growth 
rate of the ordinary hydrodynamic resistive drift instability, which are provided by shear 
flow. The right-hand side of Eq.( l4"3|) is small for (Vq) 1 < t < t s , and because of the weakness 
of ion Landau damping. We equate the right-hand side of Eq. (1431) to zero and obtain in the 
lowest order the modal solution, $ (k, t) = $o ex P ( — (k)t). Making use this solution in 
Eq. (j4"2l) to evaluate G e (k, t), we obtain 

(uj(]<.)-kyv de ) pLW ( ^( k )+^e N \ 

GAKt) = -i b ' T - V " }^:U (k,t), (45) 

kzV Te V 2 ^ y/2k z V Te J 

where W (z) = e~ z2 ^1 + (2z'/-y/7r) J e* 2 dtj is plasma dispersion function. Accounting for 
the small right-hand side of Eq.f l43|) . we seek in the next approximation the solution in the 
form 

$ (k, t) = $ exp (-ioj (k) t + a (k, t)) . (46) 

Assuming that all terms at the right-hand side of Eq. (j4ip are of the same order, by inserting 
Eq. (l46l) into Eq.( )43|) and neglecting derivative da (k, t) jdt in the right-hand side of Eq. (l43l) . 
the equation for a (k, t) is found 



t 

, _ . IdaCk, ti) 
dti$(k,*0 ^ ; 

to— >— oo 
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iScu (k) 



e(ti 



-2t l [l + 



T + k\ p 2 

26,; 



0. 



(47) 



in which was assumed that to — > — oo. In this equation, Sco (k) = KeSu (k) + 27 (k), where 
with assumption that v e ^> oj (k) and u e ^> k z vr e , 



Re<Mk) 



a; 3 (k) k\p 2 s 



2 ' 



and 



7 (k) 



1 + A;ip2) fc 2 %e 



^ 2 (k) fcjpg 

4e 



(48) 



(49) 



1 + k 2 ± p 2 ) kin 

is well known growth rate of the resistive drift instability^. We obtain for (V^') 1 < t < t 
the solution to Eq. ()471) in the form 

(1 + T) t 3 



$ (k, t) = $ exp 



-iu (k) t 



+ 700*- 



3(T + ^) t 2 6 , 



+ iRe 5w (k) t 



1 + 



2bi 



(50) 



in which effect of the shear flow is the same as it was obtained in the linear non-modal 
theory of the kinetic drift instability 13 - shear flow leads to the decrease the frequency and 
the growth rate of the resistive drift instability with time. 

In long-time limit, t ^> t s , arguments of Bessel function J in Eq. fj4"0"j) become large, 
k (t) k (ti) p 2 ~ t 2 /t 2 s > 1 and asymptotics (134)) simplifies Eq. (f40|) for potential <p (k, t), 

t t 



[1 + T) dh 



dip (k, t) 
dtt 



dtx 



d 

dtt 



t s 



to 



to 



+ ikyV d i J dti<p (k, t) 

ts 



to 



+ y?(k,t ) 



^ntt[ 



(51) 



t V27T 

where k\ = t~ 2 + k 2 v^. At t ^> t s , the right hand side part of Eq. llBTj) is proportional 
to small parameter t s /t ~ t s /to- Applying the methodology of the solution of Eq. (l36|) to 
Eq.f l51|) . we obtain the following solution: 



<p (k, t) = p exp 



1 t. 



. kyVdi 
l-t h 



7T k 2 z V 2 Tl t 2 s 



i+T) t y 2 Ki 

which reproduces qualitatively the same time dependence as of solution ([3 



(52) 
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CONCLUSIONS 



In this paper, we develop the non-modal kinetic theory of the hydrodynamic ion temper- 
ature gradient and resistive drift instabilities in plasma shear flow. In this theory, which is 
grounded on the shearing modes approach in the kinetic theory of plasma shear flowa^, the 
shear flow reveals as the time-dependent effect of the finite Larmor radius in the integral 
equation for the electrostatic potential. This effect is of principal importance for turbulence 
evolution in plasma shear flows. It consists in the interaction of ions undergoing cyclotron 
motion with inhomogeneous electric field of sheared modes, which due to their distortion by 
shear flows have time dependent wave number in laboratory frame. In Ref.— we obtain, that 
this effect is a source of the enhanced suppression of the kinetic drift instability by shear 
flow. In this paper we find, that in spite of their hydrodynamic nature, ion temperature 
gradient and resistive drift instabilities in shear flow at times t > (V^) -1 pass through the 
same linear non-modal kinetic processes of their evolution, as the kinetic drift instability. 
These processes reveal in the non-modal decrease with time the frequency and the growth 
rate of the unstable perturbations that display the universality in the linear description of 
the temporal evolution of the electrostatic drift instabilities of shear flows. 

It is important to note, that being qualitatively similar in their linear evolutionary theory, 
these two discussed instabilities have different nonlinear descriptions. The renormalized 
nonmodal nonlinear theory^, which accounts for the effect of turbulent scattering of ions by 
the ensemble of the sheared modes with randomly distributed initial phases, developed for 
the kinetic drift instability, is completely applicable to the resistive drift instability and gives 
completely the same result, determined by equation (68) in Ref.— , where the frequency and 
the growth rate are determined now by Eqs. (|48|) . ( 1491) . This nonlinear effect, which is absent 
in conventional gyrokinetic theory, consists in the scattering of the ion gyration angle by 
sheared perturbations and reveals in the reduction of the growth rate with time as 7 (V^t) 6 . 

The hydrodynamic ion temperature gradient instability considered here have the modal 
growth rate of the order of its frequency. Therefore, the methodology of the renormalized 
nonlinear theory, developed in Ref.—, which require the smallness of the growth rate in 
comparison with frequency, is not applicable for that instability and new nonlinear theory 
have to be developed. 
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Appendix A: TRANSFORMATION TO SHEARED COORDINATES FOR 
NON-STATIONARY SHEAR FLOW 

The procedure of the transformation of the Vlasov-Poisson system to the sheared coor- 
dinates, developed in Ref.— , is generalized easily on the case of the spatially homogeneous, 
but time dependent velocity shear, 

V(r, t ) = V (M) = = = Vi^e,. (Al) 

where Vq is a parameter with dimension of the velocity shear and a (t) is a function with 
dimension of time. The transition in the Vlasov equation from coordinates x, y, z and 
velocity v to convected with flow velocity coordinates v x = v x , v y = v y + V^da/ dt)x, 
v z = v z , and sheared coordinates, determined by x = x, y = y + V^a (t) x, z — z (it is 
assumed that inhomogeneous electric field, and respectively shear flow originate at time 
t = and time ordering u c {T 3> 1 is adapted, where T is characteristic time of the flow 
velocity variations) removes from the Vlasov equation spatial inhomogeneity introduced by 
shear flow. With leading center coordinates X, Y, determined in convective set of reference, 

X = x+ sin0, (A2) 

Y = y — cos <p — Vq'o (t) (X — x) , Z\ = z — v z t, (A3) 

and velocity space coordinates 

v x = v ± cos 0, Vy = y/rjv ± sin 0, = X - uj c /i (t) , v z = v z , (A4) 



where 



+ da(p_ = r* ^^ dtu (A5) 

u c dt J 
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we obtain Vlasov equation in form (|T6|) with time dependent coefficient 77. The electrostatic 
potential ip (x, y, z, t) in Eq. ffTB"]) is determined by f|T7|) . where now 

k\ {t) = (k x - V^a (t) k y f + -±-k 2 y , (A6) 

and tan 6 = k y j ^Jrj [t)[k x — V^a (t) k y ). Therefore, the results for the time dependent shear 
( 1A1I) are reproduced easily from ones obtained for the stationary shear by the replacements 
V^t on V 'a (t). For real experimental conditions Vq <C Ud\ so rj — 1 and fi = 1 and in such 
case Eq. ffToT) still has not explicit time dependence. As a sample of such analysis, we present 
here the extension of the solution Eg. (1501 . obtained with assumption of r\ = 1, on the time 
dependent velocity shear 



$ (k, t) = $ exp 



a 2 (t) A , T + k\p\ 



+ 7(k)t- - , rHlw? 1 + 



(A7) 



The performed analysis displays, that the effects of shear flows becomes appreciable when the 
time-dependent part of the wave number flA6j) . becomes dominant, i.e. when | V^a (t) | ^> 1. 
That condition may occur for the steady or for the growing with time velocity shear. For the 
velocity shear ( lAlj) oscillating with time as da (t) jdt ~ sina^, we have \V^a (t) \ < Vq/cuo 
and condition |V^'a (t) \ ^> 1 may be attained only for low frequency oscillation of the velocity 
shear, when the condition \Vq \ ^> u is met; effect of shear flow is negligible for u > \Vq\. 
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